Synchronization of Rossler Oscillators on Scale-free Topologies 
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We study the synchronization of Rossler oscillators as prototype of chaotic systems, when they 
are coupled on scale-free complex networks. We find that the underlying topology crucially affects 
the global synchronization properties. Especially, we show that the existence of loops facilitates the 
synchronizability of the system, whereas Rossler oscillators do not synchronize on tree- like topologies 
beyond a certain size. By considering Cayley trees, modified by various shortcuts, we find that also 
the distribution of shortest path lengths between two oscillators plays an important role for the 
global synchronization. 
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I. INTRODUCTION 

Synchronization is an ubiquitous phenomenon in na- 
ture, ranging from flashing fireflies in the Australian for- 
est 0, crickets chirping in unison ^ in natural systems, 
tremor in Parkinson's disease or epilepsy in medical ap- 
plications 0, laser arrays or Josephson junctions in 
physics 13, electrochemical oscillators in chemistry 
and designed synchronization in robotics. In particu- 
lar, synchronization properties of limit-cycle oscillators 
were studied in a number of papers (for a review see 
0), but even systems which are individually chaotic like 
Rossler oscillators, can synchronize under certain condi- 
tions. Rossler oscillators are sometimes treated as pro- 
totype of chaotic systems. According to a conjecture of 
Calenbuhr and Mikhailov , the behavior of Rossler os- 
cillators shows some universal features. For a certain 
class of interactions and under the influence of noise, 
clusters of synchronizing oscillators form above a cer- 
tain threshold of the coupling strength between the os- 
cillators, and, for larger couplings, after an intermittent 
phase, the whole set of oscillators synchronizes. 

Rossler oscillators were studied for different interac- 
tion schemes and on different geometries llO^. In this 
paper we study the conditions for Rossler oscillators to 
synchronize on scale-free network topologies. Scale-free 
networks seem to be realized in a number of natural and 
artificial systems like genetic or proteomic networks, the 
world- wide- web and the internet. Synchronization is cer- 
tainly one of the important dynamical processes, running 
on these networks, as it is supposed to be a necessary in- 
gredient for the efficient organization and functioning of 
coupled individual units, that, after all, lead to well coor- 
dinated behavior in time. Therefore we are interested in 
the compatibility of scale-free topologies with synchro- 
nization, in particular for the case that the individual 
dynamics is chaotic. While usually the synchronization 
transition is studied as a function of the coupling strength 
or the system size, we describe here, in addition to the 
usual approach, a transition to the synchronized phase 
as a function of the topology. As we shall see, when a 



tree becomes too large in size to allow for synchroniza- 
tion, synchronization becomes possible again beyond a 
critical threshold in the coupling, when a critical number 
of loops and shortcuts is introduced into the tree, while 
it is impossible for arbitrary couplings on a tree beyond 
a certain size. 

In section I we introduce the model and define the or- 
der parameters that are used to distinguish the phases 
with and without the condensates of synchronized oscil- 
lators. In the second section we describe details of the 
simulations and summarize the results in section III. 



II. THE MODEL 

We consider a system of N Rossler oscillators, dis- 
tributed on the nodes of a scale-free network, generated 
with the growth algorithm of Barabasi and Albert [l^ 
(see below). Each individual oscillator is described by 
the following set of dynamical equations that was origi- 
nally proposed by Rossler ^T§\ as a "model of a model" 
for descri bing the trajectory of flow, satisfying the Lorenz 
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For Lu = 1, a = 0.15, b = 0.2, c = 8.5 the system is in 
the chaotic state. Among various possibilities of coupling 
these oscillators, we choose 
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where Xi, yi, Zi are averages defined as 
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and accordingly for yi and Zi, hi denotes the degree of 
node i, A^j is the adjacency matrix i.e.Ay = 1 if i and 
J are connected and otherwise, that is the only place 
where the topology of the network enters. For ki = N—1 
and Aij = 1 for all i,j € 1, ...,N the system corresponds 
to a globally coupledpopulation of Rossler oscillators as 
it was considered in Isf . In our description the population 
is partially coupled rather than globally. It is coupled 
along the links of the scale-free network, therefore the 
driving force towards the common synchronized state is 
produced by nearest-neighbors, whose number varies ac- 
cording to the scale-free degree-distribution. Since our 
averages are still node-dependent, the stability analysis 

of 0, derived for x — ■^J2f=o^j^ iv^^ ahke,) does not 
immediately apply. For this case of global coupling with 
driving force that tries to reduce the difference from the 
common synchronized state (x, y, z), one expects a glob- 
ally synchronized stable state for e = 1, a < 1, so that all 
deviations from global averages exponentially decrease 
with time In our scheme the force drives to node- 
dependent average values over nearest neighbors whose 
number is neither regular nor N — 1, i.e. all-to-all. Nev- 
ertheless we find a result quite similar to the all-to-all 
case: a global attractor to a synchronized state as long 
as e < 1.25. The stability is evident on the level of nu- 
merical simulations. 

In order to check how the results depend on the non- 
linear terms of our coupling scheme, we made also some 
tests for the linear vector coupling defined according to 

Xi yi Zi -\- € (^Xi 

yi ^ Xt + ayi + e{yi - yi) 

Zi = b + {xi - c)zi + e{zi - Zi) , (4) 

i — 1, ...N, as it was used in [lll |. 

A. Choice of order parameters 

As a first indicator for a partially or fully synchronized 
state, we measure the histogram of instantaneous pair 
distances dij (t) between all pairs of nodes as a function 
of (simulation) time, defined by 

dij = [{xi - Xj)'^ + {yi - yj)'^ + {zi - z^)^] , (5) 

«,j = 1,...A^. A fully synchronized state shows up as 
a sharp peak in the distribution of dij, since the pair 
distances between any two nodes approach zero. No syn- 
chronization or desynchronization in the opposite case, 
lead to a broad distribution. As order parameters in the 
usual sense (varying between and 1), for the desyn- 
chronized phase and 1 for the fully synchronized phase, 
we choose two order parameters r and s, as proposed in 
[Tl|. defined in the following way. 

'^'^ = N(N-l) ^ ^i5-d^m, (6) 



and 

N N 

where <d{x) is the Heavyside function, i.e. 8(x) = 1 if 
a; > and Q{x) — otherwise. The parameter ^ is a 
small number to account for the finite numerical accu- 
racy, e.g. 5 — 0.0001, so that two states in phase space 
lying inside a sphere of radius 5 are considered as mu- 
tually being synchronized. The parameter r[t) gives the 
fraction of pairs of elements (j, j) which are synchronized 
at time t (i.e., dij < S). This fraction is one if all possible 
pairs are synchronized and zero if no pair is synchronized, 
intermediate values < r < 1 reflect partial synchroniza- 
tion. The second order parameter s{t) is more sensitive 
to partial synchronization. The second term on the r.h.s. 
of EqCI only contributes to the fraction if node i has 
no other node within a distance of S. Therefore s is al- 
ready 1 when the total number of states is partitionized in 
synchronized pairs without synchronization between the 
pairs. In general we have r < s < 1 (as it is confirmed 
in the figures below) if some elements form clusters while 
others are still isolated. From a simultaneous measure- 
ment of r and s it is possible to obtain some information 
about the partial synchronization that is usually a pre- 
cursor to the fully synchronized state. We measured in 
general all three functions dij, r, and s as a function of 
the number of iterations. 



III. MEASUREMENTS AND RESULTS 

A. Generating the topology 

For the scale-free topology we used the growth algo- 
rithm of Barabasi and Albert later referred to as 
the BA model. In each step, one node with m edges 
is added to the network. It is connected to m of the 
formerly generated nodes according to preferential at- 
tachment. In our simulations we chose m between 1 and 
10. For testing the role of the loops we used the regu- 
lar topology of a Cayley tree with z edges at each node, 
2 = 3, 6. The tree structure was then modified in var- 
ious ways as we shall see below. We also made some 
runs on a small-world topology, starting from a regular 
ring topology with k — 2 neighbors and randomly adding 
shortcuts to each node with probability p — 0.01 accord- 
ing to the algorithm proposed by Newman and Watts 



B. Choice of parameters 

For the parameters of the individual Rossler oscillators 
we chose a = 0.15, b — 0.20, c — 8.5 throughout all 
simulations to make sure that the individual systems are 
in the chaotic regime. The total number N of oscillators 
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FIG. 1: Histogram of all pair distances between oscillators on piG. 2: Order parameters (a) r and (b) s as a function of the 
a BA network with (a) m = 1, (b) m = 2 for various e and coupling strength e for different values of m, iV = 200 
iV = 200 



was varied between 10, 50, 200 up to 500 on the scale- 
free topology, and N — 190 on the Cayley tree. The 
parameter m of the growth algorithm varied between 1 
and 10, the coupling strength e was out of the interval 
[0.1, 1.2]. In the numerical simulations of Eq|21 we used 
the fourth order Runge-Kutta methods with a typical 
time-step size of dt — 0.001 (when N = 200). Variation 
of dt between 10^^^ < dt < 10~^ led to qualitatively the 
same results. 



C. Results for m=l 

Fig.l a) displays the results for the histogram of dis- 
tances for m = 1, iV = 200 and various couplings e up to 
1.25, above which the numerical integration becomes un- 
stable. The distributions are broad and do not indicate 
any synchronized state. This result is further supported 
by measurements of r and s as shown in Fig. 2 a) and 
b), respectively, r stays zero for m = 1, while s increases 
from e = 0.1 on, indicating some partial synchronization. 
The value of TV = 200 seems to represent the large- 
limit, for the considered range of e and to, since we ob- 
tained the same result for r and s for TV = 300, 400, 500. 
On the other hand, for smaller systems, iV < 20, we 
do see a fully synchronized state when the coupling e 
exceeds a critical threshold. As value for i5 in Ea. (|6l7|l 
that accounts for the finite numerical accuracy, we choose 



6 = 0.0001. For too small values of 6 we observe large 
variations in the long-time behavior of r and s, for too 
large 6, the values of r and s are stable over long times, 
but their values do depend on S. In between, i.e. for 
0.0001 < S < 0.1 we find a plateau for the values of r 
and s, that is, the behavior becomes independent of the 
size of 6. 



D. Results for m larger than one but still integer 

If we keep the number N of oscillators fixed to 200, 
we observe for to > 1 a fully synchronized state above a 
critical threshold in the coupling e; this threshold is the 
larger the smaller to, again r < s in general, as seen from 
Fig. 2 a) and b). 



E. Results for intermediate noninteger m 

One of the main differences between the Barabasi- 
Albcrt networks with parameter to = 1 and m > 1 is 
the tree-like structure for m = 1 and the existence of 
loops for TO > 1. In order to check whether it is only 
the loops that facilitate synchronization and how many 
loops are needed, we generalized the growth algorithm to 
non-integer values of m in the following way. We intro- 
duce an additional probability pm for a new node to have 
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FIG. 3: Histogram of pair distances on a BA-model with 1 < 
(m) < 2 with N = 200 and e = 0.9 




FIG. 4: Order parameters r (a) and s (b) as function of the 
parameter m that is used to distinguish different topologies, 
for two values of the coupling strength 



TO = 1 edges and 1 —pm for having to = 2 edges attached 
to the nodes of the network when it is introduced during 
the growth process. The distribution of pair distances 
of oscillators for 1 < (to) < 2 is displayed in Fig|31 For 
given N and e we therefore observe a transition to a fully 
synchronized state as a function of " topology" , parame- 
terized via the parameter to. For N = 200 and e — 0.9, 
we have 1.35 < TOc < 1.4, cf. Fig|31 Figgja) shows that 




FIG. 5: Cayley tree for 2 = 3 and additional interconnections 
(dashed lines) attached according to the three rules a), b), c), 
respectively, as described in the text. 



the position of the transition, now in to rather than in 
e, depends on the coupling strength for fixed N. The 
smaller e, the larger TOc. An interesting feature is seen in 
Fig^Ja), where s is plotted as a function of to. For e = 0.3 
and m ~ I, the finite value of s indicates some partial 
synchronization, s then drops to zero at to = 2 and in- 
creases to 1 for TO = 3. As we have argued above, s = 1 
does not necessarily imply full synchronization, but some 
partial one, at least. The behavior of s is non-monotonic 
as a function of to. A similar non-monotonic behavior of 
s as function of time was observed in "o] for Rossler oscil- 
lators, for which a partial synchronization was followed 
by desynchronization, before the full synchronization set 



F. Rossler oscillators on a Cayley tree 

From the former results we conclude that a certain 
number of loops facilitates synchronization on scale-free 
networks, the larger to, the more loops 16], the smaller 
the coupling strength needed for synchronization. In or- 
der to check whether it is the mere number of loops that 
facilitate synchronization or also the type of loops, we 
studied N Rossler oscillators on a Cayley tree whose 
regular structure was modified in a controlled way by 
adding a) edges to construct a given number of triangles 
at random locations, b) shortcuts between the outermost 
nodes, c) shortcuts between the outermost nodes and the 
central node with probability pi, see Fig|51 For > 187 
and arbitrary values of e (more precisely, for e as large 
as e = 0.9), Rossler oscillators do not synchronize on 
a Cayley tree, neither for shortcuts according to a) or 
b), but for method c) and pi > 0.9 they do so, as it is 
seen in Fig|H| As it turns out also from measurements 
on a small- world topology with p — 0.01, synchroniza- 
tion in all of our measurements goes along with a dis- 
tribution of shortest path lengths that looks Poissonian 
like, cf. FigEt>) On the other hand, it is not the average 
shortest path length that is conclusive alone. As Fig. [7| 
shows, the average is almost the same in a) and b) while 
the distribution is different and those of Fig[7Ji) corre- 
spond to desynchronization on a Cayley tree without or 
with shortcuts between the outermost nodes (apart from 
Pi = 0.99), while those of Fig|7|D) lead to synchronization 
apart from the BA model with m — I, since in this case 
the loops are totally absent. This shows that neither the 
existence of loops alone nor the existence of a certain dis- 
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tribution of shortest path lengths alone are sufficient to 
guarantee synchronization. 



IV. SUMMARY AND CONCLUSIONS 



FIG. 6: Distance distributions for Rossler oscillators on a 
Cayley tree with additional edges according to c) with differ- 
ent interconnection probabilities pi, as further explained in 
the text 
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FIG. 7: Histograms of shortest paths of length £ for vari- 
ous networks with Rossler oscillators. The envelope of the 
distributions looks similar for a) with p = 0.99 and b) for 
m = 1, m = 2 and on the small world network 



We studied an ensemble of Rossler oscillators on sale- 
free networks constructed by the Barabasi- Albert growth 
algorithm. In contrast to the usual investigations we 
studied the transition from the desynchronized or par- 
tially synchronized state to the fully synchronized state 
as a function of the network topology, parameterized by 
the m, the number of newly attached edges in the growth 
algorithm. For the tree topology (m = 1) and given cou- 
pling strength e, there is a fully synchronized state below 
some critical size A'^ that disappears for larger N. This re- 
sult is similar to synchronization of Kuramoto oscillators 
on Cayley trees which is possible for small enough size N 
and coordination number z |l7j . Above a certain number 
of nodes, the tree of Rossler oscillators is no longer syn- 
chronizable however large the coupling strength is, but 
it is then the parameter m that introduces loops and 
shortcuts into the tree, and along with this allows for full 
synchronization again when m exceeds a certain value 
that depends on A^ and e. The threshold in e depends on 
A^ and m, and vice versa, the threshold in m is sensitive 
to A^ and e. Small A^, large e (chosen from the stability 
regime) and large m favor synchronization. These quali- 
tative results are not specific for our choice of nonlinear 
couplings between the Rossler oscillators, but also hold 
for the vector coupling scheme of Eq2| Moreover, numer- 
ical simulations of Rossler oscillators on Cayley trees with 
difi^erent artificially introduced shortcuts suggest that it 
is not only the mere number of loops that favors syn- 
chronization, but loops that provide real shortcuts such 
as those between outermost and central nodes. 
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